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' Abstract 

| We consider the effect of a random magnetic field in the convective zone of the Sun 

^C) ■ superimposed to a regular magnetic field on resonant neutrino spin- flavour oscillations. 

0^ ! We argue for the existence of a field of strongly chaotic nature at the bottom of the 

convective zone. In contrast to previous attempts we employ a model motivated 
regular magnetic field profile: it is a static field, solution to the solar equilibrium 
hydromagnetic equations. These solutions has been known for a long time in the 
literature, we show for the first time that in addition they are twisting solutions. In this 
scenario electron antineutrinos are produced through cascades like v e L — > — > v e R, 
O I The detection of v e R at Earth would be a long-awaited signature of the Majorana 

nature of neutrino. The expected signals in the different experiments (SK,GALLEX- 
SAGE,Homestake) are obtained as a function of the level of noise, regular magnetic 
field and neutrino mixing parameters. Previous results obtained for small mixing and 
ad-hoc regular magnetic profiles are reobtained. We confirm the strong suppression 
for large part of the parameter space of the f^-flux for high energy boron neutrinos 
d ' in agreement with present data of the SK experiment. 

We find that MSW regions (Am 2 10~ 5 eV 2 , both small and large mixing solu- 
tions) are stable up to very large levels of noise (P=0.7-0.8) but they are acceptable 
from the point of view of antineutrino production only for moderate levels of noise 
(P« 0.95). 

For strong noise and reasonable regular magnetic field, any parameter region 
(Am 2 , sin 2 29) is excluded. As a consequence, we are allowed to reverse the prob- 
lem and to put limits on r.m.s field strength and transition magnetic moments by 
demanding a particle physics solution to the SNP under this scenario. 
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1 Introduction 



A neutrino transition magnetic moment can account, both for the observed deficiency of 
the solar neutrino flux and the time variations of the signal. The overall deficit is caused 
by the alteration or suppression of the neutrino energy spectrum. The time dependence 
may be caused by time variations of the magnetic field in the convective zone of the Sun. 
As it has been shown in |lj, magnetic moments solutions are not ruled out by present 
experimental data; that is valid for both, absolute deficits and time variations of the 
observed solar neutrino flux. 

In transverse magnetic fields, neutrinos with transition magnetic moments will expe- 
rience spin and flavour rotation simultaneously (resonant spin-flavour precession, RSFP). 
The observation of electron antineutrinos from the Sun would lead to the conclusion 
that the neutrinos are Majorana particles. There are however stringent bounds on the 
presence of solar electron antineutrinos coming from the High energy Boron neutrinos ( 
E v >«8.3 MeV g §). 

Magnetic fields measured on the surface of the Sun are weaker than within the interior 
of the convective zone. The mean field value over the solar disc is of the order of 1 G and in 
the solar spots magnetic field strength reaches 1 KG. In magnetic hydrodynamics (MHD) 
one can explain such fields in a self-consistent way if these fields are generated by dynamo 
mechanism at the bottom of the convective zone. In this region the strength of small scale 
regular magnetic fields could reach a value of 100 KG. These fields propagate through the 
convective zone and photosphere decreasing in the strength value while increasing in the 
scale giving traces in from of loops in bipolar active regions (solar spots). 

Large-scale toroidal magnetic field created by dynamo mechanism in convective zone 
has strength even less than small-scale r.m.s fields near the bottom of convective zone. This 
is the main reason why one should consider neutrino propagation in the random magnetic 
field of the Sun. The ratio of the r.m.s. random field and the regular (toroidal) field may 
be ~ 20-50, therefore the problem of RSFP neutrino propagation in noisy magnetic field 
seems to be important. 

Estimations for the ratio of rms fields to regular field are necessarily very rough. In 
textbooks [|], ||, ^, |7j] we find the conservative ratio 

(&)/Bl ~ 1. 

In more elaborate models the ratio of magnetic energy densities is given by the magnetic 
Reynolds number, 

{B 2 )/Bl ~ R^, 

which may be much bigger than unity for plasma with large conductivity. Here a > is 
a topology index [|, |, |, §. 

The effect of random magnetic fields in RSFP solutions to the SNP and antineutrino 
production has been explored previously for simplified models [H]. In this work we 
will deal with the complete problem, we will present calculations of neutrino spin flavour 
conversions in presence of matter and magnetic field. The magnetic field will have two 
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ingredients. The first ingredient will be a a theoretically motivated solar magnetic field 
profile. This magnetic field, which is one solution to the static magnetic hydrodynamic 
equations, is in addition a twisting field. The degree of twisting, the transversal profile 
and the ratio core field/ convection field are functionally related under this model. Of 
all possible solutions, in practical calculations we have considered only the solution with 
minimum twist which is at the same time the solution which implies a minimum value for 
the central field ( a factor 2 or less of the magnetic field at the convective zone). As a 
second ingredient, the effect of a layer of magnetic noise generated at the bottom of the 
convective zone has been included, we have justified that the level of noise in this region 
can be certainly very high. 

Notice that an analogous problem but without the regular twist field motivated here 
while with use of a direct numerical method for dealing with the random magnetic fields 
in the solar convective zone has been recently considered ([12|). The results obtained in 
fl2|| are in general consistent with those obtained in the present work (see Conclusions 
below) . 

In Section (2) we present the effective Hamiltonian governing the evolution of the 
4x4 neutrino wave function (two flavors times two helicities) in presence of a generic 
magnetic field. In Section (3.1) we give some physical arguments supporting the idea 
that r.m.s fields could be indeed important in the solar convective region. The analytical 
expressions for the hydro dynamic-inspired regular large-scale magnetic field appear in 
Section (3.2). There we show how different configurations of this field can be classified 
according to its intrinsic "twist" in the perpendicular plane to neutrino propagation. In 
Section (4) we give the master equation governing the time evolution of the averaged 
density matrix in presence of the random magnetic field. We show analytically how in 
the limit where chaotic field dominates, the 4x4 evolution equations decouples and the 
averaged transition probabilities follow a Markovian process. The results of the numerical 
integration of the exact averaged master equation are shown and discussed in Sections (5) 
and (6). Prom the calculated transition probabilities the expected total signals in each 
of the existing experiments (SK, SAGE-GALLEX,HOMESTAKE) are obtained and x 2 - 
allowed parameter regions. In addition the expected electron antineutrino signal in SK is 
calculated and compared with existing bounds. 



2 The Master equation 

There are two channels for v e L — ► f e _R-conversions corresponding to the cascades: (i) 
v e L V[iL v e R or (ii) u e L — > v^n — > D e R. The former case realizes if we assume 
a zero field in the radiative zone and in the core or if we exclude right-handed neutrino 
production there and the MSW-conversion [y e t — > v^l) takes place before neutrino reaches 
convective zone. The latter case (ii) realizes when a strong magnetic field is present both 
in the radiative zone and in the solar core. In this case the RSFP takes place before the 
MSW-conversion since the v e L — v^r energy splitting V = Gi?y / 2(A r e (ri) — N n (n)) is less 
than for the MSW one (u b l — ► v^i) at the same point V = GpV2N e (ri). This is true for 
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the typical admixture of helium and heavy elements considered to be present in the inner 
layers of the Sun. 

We consider conversions v e L.R - > v aL.R-, a = p or t, (for definiteness we will refer to p 
for the rest of this work) for two neutrino flavors obeying the master evolution equation 

/V e -c 2 6 s 2 5 pB+{t) \ 

-l*BZ(t) 



id, 



VeR 
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v^l s 2 5 -fiB'Kt) V^ + c 2 5 v^l 

\VpB.) V/x^lO) s 2 5 -V^ + c 2 5/ \V pR J 

where c 2 = cos 29, s 2 = sin(2#), 5 = Am 2 /4E are the neutrino mixing parameters; p 
is the neutrino active- active transition magnetic moment; 

Bf(t) = B± ± (t)+Bf(t) (2) 

is the magnetic field component which is perpendicular to the neutrino trajectory in the 
Sun; 

V e {t) = G F V2(p(t)/m p ){Y e - Y n /2) 

and 

V„(t) = G F V2(p(t)/m p )(-Y n /2) 

are the neutrino vector potentials for v e L and v^l in the Sun given by the abundances of 
the electron (Y e = m p N e (t) / p(t)) and neutron (Y n = m p N n (t) / p(t)) components and by 
the SSM density profile @ 

p{t) = 250 gem- 3 exp(-10.54 t). (3) 

The transverse magnetic field B±(t) appearing in Eq.(Q) is given by the following 
expression where Cartesian and polar coordinates are written explicitly: 

Bf(t) = B x (t) ± iB y (t) = \B ± (t)\e± m ® (4) 

where 

= arctan B y / B x . (5) 

The Cartesian form is useful in writing the averaged master equation. The polar form is 
the most convenient in separating out the twisting part. 

The nature and magnitude for the regular (-BojJ and chaotic parts (B±) of the magnetic 
field will be the subject of the next section. The equation for the evolution of the average 
density matrix corresponding to the master equation (|l]) can be found using the formalism 
developed in and will presented later in this this work. 

The magnetic field strength enters the evolution Eq.(^) being multiplied by the neutrino 
transition moment p. The existing upper limits on the magnetic moment of the electron 
neutrino include the laboratory bound p < 3 — 4 x 10~ 10 /i£ from reactor experiments as 
well as stronger (one or two orders or magnitude) astrophysical and cosmological limits. In 
our calculations we will consider always the product pB. Expected values of B ~ 1 — 100 
kG in the Sun convective zone and p = 10~ n /i£ would give an expected range for the 
product pB ~ 10~ 8 — 10~ 6 psG ~ 5.6 10~ 17 — 10 -15 eV or in the practical units which 
will be used throughout this work pB ~ 0.1 — 10.0 puB^. 
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3 Solar magnetic fields 



3.1 Random magnetic fields 



The r.m.s. random component y (B 2 (t)} can be comparable with the regular one, Bo(t), 
and maybe even much stronger than Bq, if a large magnetic Reynolds number R m leads 
to the effective dynamo enhancement of small-scale (random) magnetic fields. 

Let us give simple estimates of the magnetic Reynolds number R m = lv/u m in the 
convective zone for fully ionized hydrogen plasma (T 3> Ih ~ 13.6 eV ~ 10 5 K). Here 
I ~ 10 8 cm is the size of eddy (of the order of magnitude of a granule size) with the 
turbulent velocity inside of it v ~ va ~ 10 5 cm/s where va = Bo/y/Anp is the Alfven 
velocity for MHD plasma, Bq is a large-scale field in convective zone and p is the matter 
density (in g/cm 3 ) in the SSM. 

The magnetic diffusion coefficient Vm — 

c r /AircTamd (~ magnetic viscosity) enters the 

diffusion term of the Faraday equation, 

^> = rot[v x B(t)] + u m AB . (6) 
at 

Here c is the light velocity; the conductivity of the hydrogen plasma cr C ond = ^pi/^^ep- 



Upi = y4-7re 2 n e /m e = 5.65 x 10 ^Jn~ e s 

is the plasma (Langmuir) frequency; u ep = 50n e /T 3 / 2 s~ 1 is the electron-proton collision 
frequency, the electron density n e (= n p ) (cm -3 ) and the temperature T (K). 
Thus we find that the magnetic diffusion coefficient 

10 13 (T/1 K)-^ 2 cmV 1 



■'m 



does not depend on the charge density n e and it is very small (the Reynolds number is 
big) for hot plasma T > 10 5 K S> 1 K. Actually, from the comparison of the first and 
second terms in the r.h.s. of the Faraday equation Eq. (^) we find that v/l 3> v m /l 2 , or 
v m <C vl ~ 10 13 cm 2 s~ 1 since T/l K S> 1. This means that the magnetic field in the 
Sun is mainly frozen-in. Neglecting the second term in Eq. (^) and using the Maxwell 
equation 

rotE = - C - l {dB/dt) 
we obtain the condition for frozen-in field: the Lorentz force vanishes, 

~ (E+ [vx B]/c) « 

but the current 

J = cr cond (E + [v x B]/c) 
remains finite if the conductivity is large, a C ond oo- 
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The magnetic Reynolds number 

R m = lvuj 2 pl /{c 2 v ep ) ~lvx 10" 13 (T/1 Kfl 2 cnJs' 1 

is huge if we substitute the estimate Iv ~ 10 13 cmPs^ 1 given above. A large value for 
the Reynolds number is a necessary condition for an effective dynamo enhancement in the 
convective zone. 

The estimation of the quantity n for the solar convective zone (and other cosmic 
dynamos) is the matter of current scientific discussions. The most conservative estimate, 
simply based on equipartition, is n = constant. According to direct observations of galactic 



magnetic field presumably driven by a dynamo, n ~ 1.8 [15]. A more developed theory of 
equipartition gives, say, n 47rlni? m (sec fl6|]). Notice that this estimate is considered 
now as very conservative. Basing on more detailed theories of MHD turbulence estimates 
like b ~ yJR m B are discussed || , |J . 

The random magnetic field component in the Sun {{B{t)) = 0) will be described in 
general by an arbitrary correlator 

(B(t)B(t)) = (B 2 )f(t-t). 

We will assume that the strength of the r.m.s. field squared (B 2 ) = i]Bq is parametrized 
by the dimensionless parameter n = > 1 , which it can be, in general, much bigger 
than unity, n 3> 1. 

The correlator function f(t) is unknown a priori but it takes the particular ^-correlator 
form f(t) = Lo5(t) if the correlation length (for two neighboring magnetic field domains) 
is much less than the neutrino oscillation length, Lq <C lose- Lq can be considered a 
free parameter ranging in the interval 1 — 10 4 km. In the averaged evolution equations 
it appears only the product t)Lq. Thus in what follows we will present our results as a 
function of the quantity P which is a simple function of such product: 

P = ~(l + exp(- 7 )) (7) 

7 = |o 2 Ai = ^nL (uB ) 2 At. (8) 

The reason for using P is that it is a good approximation for the depolarization that 
the presence of noise induces in the averaged neutrino density matrix. At is the distance 
over which the noise is acting. We have supposed in our computations that the noise is 
effective only in a thin layer with thickness At = 0.1 R@ starting at r = 0.7 Rq, the 
bottom of the convective zone. This is represented together with the regular transverse 



profile in Fig. ([T]) . In Table @ the quantities y (B 2 ) and rj are computed for a given P 
supposing the reasonable value Lq = 1000 Km. 

3.2 Regular large-scale magnetic field in the Sun. Twist field. 

Many phenomenological formulas for Bq± in the convective zone and in the central region 
of the Sun p?]], with or without twist, has been employed in the literature. In particular, 
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twist was applied for effective enhancement of the process v e L — > z> e # in the Sun (i.e. in 
Ref . fig] ) . Without twist the coherent sum of two amplitudes for the process v e i — > v e R, 



Mi(v eL -> f Mi -> z> ei? ) + M 2 (v eL -> i^ij -> i> 6jR ), 

occurs proportional to the small abundance of neutrons and production of electron an- 
tineutrino is negligible [Qj. 

In the presence of the twist, <$(i) 7^ 0, this sum becomes proportional to the angular 
velocity 

$ = kR- 1 = K x 0.3 x 10~ 15 eV 

where [k] = n = 1, 2, ... is the effective number of revolutions of the field -Bj_ in the plane 
which is perpendicular to the neutrino trajectory. The more rounds (n>l) take place 
before the resonant RSFP position, given by the condition 

G F \/2(N e - N n ) - 25 cos 2(9 -6 = 0, 

the deeper the resonant point happens (for successive (positive) sign of <3?). It is possible 
then the merging of RSFP [v^l — ► u e R) and MSW (u e L — > v^£) resonances since the MSW 
resonance is unchanged, {Gp^/2N e = 25 cos 29). For the case of merging of resonances the 
of v e L — ► ^eij-conversions tend to be more adiabatical and complete than for separated 
resonances |I~9| ], 

Unfortunately, in the known MHD plasma solutions for the toroidal field evolving in 
convective zone due to dynamo mechanism (Yoishimura model, [20]) the predicted number 



of revolutions for twist (around toroids in northern and southern hemispheres of the Sun) 
is very small: n > 1. For this small twist rate $ ~ n the resonances above are fulfilled 
for a small Am 2 ~ 10 -7 — 10 -8 eV 2 and happen not too deep in solar interior to have 
important consequences. 

Below we apply for the neutrino conversions described by our master equation Eq. 
(|l|) the self-consistent model of large-scale regular field given in |ffif |. The global solar 
magnetic field is the axisymmetric equilibrium solution of the MHD static equations (quiet 
Sun) in the spherically symmetric gravitational field of the Sun. The reasonable boundary 
condition Bq = on the photosphere ( r = Rq) is imposed in addition. Any field solution 
to these equations and boundary conditions is a twisting field with an arbitrary small or 
large number of revolutions along radius (k = 1, oo, the twist rate can be taken as label 
for distinguish particular solution within the family). 

The spherical components of the magnetic field for the fundamental radial mode n = I 
do not depend on azimuthal angle (p. In the whole region < r < R & they are of the form 
(the different solutions are labeled by k): 



B 0r {r,9 s ) = 2Kcos9 s l-- 2 ^h^i 

L r^Znu sin Zziu V 



Boe a (r,Os) = -Ksm6 s 



r 2 z 2 k sin2z 2 k v ar 
3R%, /sin ar 



cos ar 



3-Rq / sin ar 
a 

2 H — ^ (- — — cos ar — ar sin ar 



r 2 z 2 ksm2z 2 k^ ar 
r 3-Rq / sin ar 



B <f>(r,9 s ) = Kz 2k sm9 s — r^z cosar . (9) 

LRq rz 2k sm2z 2k \ ar > J 
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Here 9 S is the polar angle. z 2 k is any root of the spherical Bessel function 

fn(z) = \fzJ n+1 / 2 {z) 

of the first order (n =1). This equation follows from the boundary condition. For the first 
three roots: 

z 2k aR Q = 5.7635, 9.0950, 12.3224,... (k = 1, 2, 3, ...). 
The constant K is related with central field B core , the only free parameter in this model: 

K = - ^ -. (10) 

2(1 -<xR / sin aRo) 1 ; 

The modulus of the perpendicular component is of the form: 



B ± = jB*+Bl e =B co J—f(r) (11) 



r 



where f(r) is some known function of gentle behavior. Obviously, on the solar equator 
(9 S = 0) the perpendicular component vanishes. This is the same behavior as it is shown 
by the toroidal field in the Yoishimura model [£(]] El- 
Notice the regular behavior of Bq at the Sun center (r = 0): 

B r (0) = B core cos9 s , (12) 
Be s (0) = -Bc^-e sin S , (13) 
B^O) = B core sm6 s r/R Q ^0. (14) 

The twist for the perpendicular component B$^{r) is defined by the angle: 

®(r)z2k = arctan B^r) /B (r) ~ z 2 k (15) 
6(r) 22fc « kR^. (16) 

For k = 1 for example, changes between and ir along the neutrino path from the 
core of the Sun through the surface at an uniform rate. For higher roots we observe that 
effectively $ ~ kn uniformly. 

According to this model, the expected magnetic field at the core is typically only 2- 
3 times (or less) the magnetic field at the convective zone. For the values that we will 
consider later, l?o.7 ~< 100 — 200 kG, the values corresponding at the core are well below 
astrophysical bounds derived from traces of these fields at solar surface. 



x This property is general and could lead to semiannual variations of the neutrino flux if one takes into 
account inclination (7°) of the solar axis (perpendicular to solar equator) to the ecliptic with the Earth 
orbit. 
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4 The averaged master equation. 

The master Equation ([!]) can be written in terms of the density matrix p(t) as: 

id tP = [H reg ,p\+pB x {t)[V x ,p]+pB y (t)[V y ,p\. (17) 

The elements of the matrices Hq, V x , V y can be read off the Eq. (|]). The B x , B y are the 
Cartesian transversal components of the chaotic magnetic field. Vacuum mixing terms 
and matter terms corresponding to the SSM density profile given before and the regular 
magnetic part Hamiltonian, with the profile determined by the Eqs.(^), are all included 
in H reg . In particular, the matrices V x>y are given in terms of the Pauli matrices a\p, by: 



It is our objective in this section to write the differential evolution equation for the 
average density matrix (p) . 

We assume that the components B x , B y are statically independent, each of them char- 
acterized by a delta-correlation function: 

(B x>y (t)B x , y (t')) = (Bl y )L 5(t-t) (19) 
(B x (t)B y (t')) = 0. (20) 

From now on we will make the following assumption, which is reasonable from equipartition 
arguments, 

(B 2 X ) = (B 2 ) = (B 2 01 _)/2 = (B 2 )/3 
The averaged evolution equation is a simple generalization (see Ref. ]l4j for a complete 



derivation) of the well known Redfield equation 23] for two independent sources of 
noise and reads (ft 2 = L p 2 (B 2 }/2 = r]L Q (pB ) 2 / '3): 

id t (p) = [H reg , ( P )] - in 2 [v x , [v x , ( P )]] - in 2 [v y , [v y , ( P )]]. (21) 



It is possible to write the Eq.(|2lD in a more evolved form. Taking into account the 
particular form of the matrices V X)V we can simplify considerably the double commutators: 

id t (p) = [H reg , (p)] - 2in 2 (2(p) - V x (p)V x - V y (p)V y ) . (22) 

The second term in Eq.(^), which gives the leading noise behavior, can be eliminated 
by a rescaling of the density matrix: 

(p(t))=exp(-An 2 t)(p'(t)). 

With this redefinition the evolution equation reads: 

id t (p') = [H reg , (pf)] + i2Q 2 (V x (p')V x + V y (p')V y ) . (23) 
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The expression in Eq. ( p3[ ) is the one which has been used in the numerical calculations 
to be presented below. It is useful however to consider the solution to Eq.(|23|) when 



H r eg = 0. This is the appropriate limit when dealing with extremely low Am 2 or very 
large energies, for an extreme level of noise or when the distance over which the noise is 
acting is small enough to consider the evolution given by H reg negligible. In any other 
scenario it can give at least an idea of the general behavior of the solutions to the full 
Eq.(^). When Hq = only the two last terms in the equation remain and an exact 
simple expression is obtainable by ordinary algebraic methods. The full 4x4 Hamiltonian 
decouples in 2x2 blocks. The quantities of interest, the averaged transition probabilities, 
are given by the diagonal elements of (p) . If Pf^ are the final and initial probabilities ( at 
the exit and at the entrance of the noise region) their averaged counterparts fulfill linear 
relations among them, schematically: 

Qf B = MQ?> B (24) 
with Q A ' B any of the two dimensional vectors 

Q A = ( f P f eL ~ eL l ) Q B = ( Tp( eL " DeR l ) ( 25 ) 

V (P{y eL v^r)) ) \ {P(v eL vyj,)) J 

and the Markovian matrix M: 

M = (if P 1 ~ P P ) (26) 

with P defined in Eq.(^). It can be shown that in this simple case P is exactly the final 
polarization of the density matrix. In the general case with a finite H reg it can be shown 
numerically that the quantity P still gives a reasonable approximation (< 10%) to the 
real polarization, at least for the cases of interest in this work. 



5 Results and Discussion. 

The present status of the Solar neutrino problem represented by the ratios between the 
solar neutrino observations and the solar model predictions is summarized in Table (1) of 
Ref.([^]) which for the sake of completeness we reproduce in (|lj). It is claimed in addition 
by a variety of analysis (check Ref.^5| for an useful review) that Homestake, the chlorine 
radiochemical experiment, presents an anticorrelation with the solar activity. Signal time 
variations are not observed by the Kamiokande experiment. The allowed time variation 



is restricted by to be «< 30% at 90% C.L. [26 



We have calculated the expected neutrino signals in the Homestake, Ga-Ge and (Super)- 
Kamiokande experiments. For this objective, the time averaged survival and transition 
probabilities have been obtained by numerical integration of the ensemble averaged master 
equation (|i~7|) for a regular magnetic profile Bq± as given in Eqs.(|9,ri) for k = 1 and the 
matter density given by Eq.(H). 



10 



The free parameters of our model are four: the squared mass differences (6 = Am 2 /2E), 
the flavour mixing angle (s| = sin 2 29), a noise strength parameter (P) and the product 
of magnetic field and moment (nB ± at a given radius and latitude). We have found 
convenient to use the magnitude of the transverse magnetic field at the bottom of the 
convective zone (r = 0.7 R & ) for a solar latitude 9 S = 7°; for other latitudes they have to 
be rescaled accordingly. 

In Figs.(^-^) we show the electron neutrino survival and transition probabilities for 
some illustrative cases. In Fig.(|2|) the dependence of the probabilities with Bq± for a 
fixed small mixing angle (s 2 = 0.01) and in absence of noisy magnetic field is shown. In 
general spin-flavour precession is suppressed by large mass differences, the non-suppression 
condition is \iB >~ S. One can observe the effect of this suppression around 5 ~ 10~ 7 — 
10~ 8 eV 2 /MeV. Magnetic precession and MSW transitions respectively predominate below 
and above this value. For the range of values considered in this case we observe a very 
modest production of F e 's (always smaller <ss 0.5% ). The production of F^'s can be 
however very important. In Fig.(||) we show the dependence with s 2 for a fixed magnetic 
field of very low magnitude (fiB = 1.0 fJ-nB^ = 10 _1 Vb x 10 kG). We observe again the 
same transition region around 5 10~ 7 - 10~ 8 eV 2 /MeV, below this range the magnetic 
precession dominates but for this low value of the magnetic field the conversion probability 
is anyway very small. No significant quantity of either electron or muon antineutrinos is 
produced irrespective of the mixing angle. For stronger magnetic fields (of the order 
fiB ~ 15 /U11-B4) the results are similar, the triangular regions, distinctive structure of the 
MSW effect, appear strongly distorted, the production of muon antineutrinos is notable 
for 5 ~ 10~ 6 — 10 -9 eV 2 /MeV with mixing angles s 2 . < 0.6 — 0.7. The production of 
electron antineutrinos is not significant even at these high magnetic fields. 

In the next two figures we show the expected signals in the Chlorine, Gallium and 
Kamiokande experiments for a variety of parameter combinations and fixing the transverse 
magnetic field to a relative high value: fiB = lbnnB^. Figs.(§,^) correspond respectively 
to the cases P=1.00 (absence of magnetic noise), P=0.70 (strong noise). The signals 
in each of the experiments have been normalized to the value in absence of neutrino 
oscillations. Comparing both figures we can observe in the noisy case how the values for 
the signal approach 1/2. This will be important later when we compute the parameter 
areas consistent with observed total rates, this smoothening, particularly in the Ga-Ge 
experiment will be the responsible for the appearance of allowed regions at high magnetic 
fields. 

For the Kamiokande experiment we have taken into account the four elastic scattering 
processes u x e — > v x e with v x = z/ e , z7 e , z7„. Using the expressions for the elastic cross 



sections appearing in [27] we have defined an "effective" conversion probability: 

P eff,e = E C *^ (27) 
x 

C-x — &x\Ev> Te ) min)/'^u e (E u , Te,min) (28) 

The coefficients c x depend on T e rnin > the kinetic energy threshold for the observation 
of the scattered electron. For T e m i n MeV the coefficients c x are very accurately 
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independent of the neutrino energy Ey for values above ~ 1 IVteV. For 2~e 5 rmn 

« 7 MeV 

the coefficients associated with the muon neutrino and antineutrino are still practically 
energy independent while c$ e is only slightly energy dependent. 

The effective conversion probability is defined in such a way that the signal expected 
at the Super-Kamiokande experiment is: 

Ssk = N J dE e(E)a{E) v ME)P ef fA E ) ( 29 ) 

where e(= 1. for us) is the experimental detection efficiency, a Ve the electron neutrino- 
neutrino elastic cross section and $ the total incoming electron neutrino flux predicted by 
SSM. N is a normalization constant. 

The weighted electron antineutrino appearance probability is defined as: 

, p v I E >EtkdE *(E)<S>(E)P Ue , Ve (E) 

{ veVe) ~ J E > Eth dE a(EME) ^ 

where a(E),<&(E) are respectively the differential cross section for the process isotropic 
background V e + p — ► e + n and the differential total neutrino flux coming from the Sun 
according to the SSM The threshold energy has been chosen as E t h = 9 MeV. From 
Kamiokande data we have the bound {P Ue v e ) =ss 0.05 ]2|, [3|, . 

In Figs.(|(||]) we present the (Am 2 , sin 2 29) exclusion plots from a combined x 2 analysis 
of the three experiments corresponding to the expected signals showed previously. 

First we comment the results in complete absence of noise (P=l) which are represented 
in Fig.@. For negligible or low regular magnetic field we observe the high squared mass 
difference solutions proportioned by the matter MSW effect. As the magnitude of the 
magnetic field increases new solutions appear and disappear in a complicated manner. 
The high angle MSW solution rests practically unmodified for all the range considered. 
The low angle solution however disappears at high magnetic fields, after experiencing 
some distortion coming from its merging with newborn magnetic solutions (compare low 
angle allowed regions in Plots (B) and (C)). The antineutrino production (dashed lines) 
is in general low and Kamiokande bounds are not specially restrictive except at very 
high magnetic fields. Note in Plot (D) the very different behavior of the two existing 
allowed regions: while the in MSW region the antineutrino production is in the 0.1 — 1%, 
compatible comfortably with Kamiokande bounds, the RSFP solution reach a value well 
above 10% and is excluded by them. It seems apparent that there are acceptable particle 
solutions to the SNP even for very large regular magnetic fields. 

The pattern of the electron antineutrino probability is very different when a small 
level of noise (P=0.95, Fig.(^)) is switched on. For Am 2 > 10~ 6 eV 2 : the antineutrino 
iso-probability lines follow the characteristic MSW triangular patterns in this region. The 
structure of the allowed regions remain unmodified. The electron antineutrino yield in 
these regions is below the 5% level for all cases except in the same region at high mixing 
angle and low squared mass difference as before. Note that only low mixing angle solu- 
tions with moderate regular magnetic fields would be acceptable if future data situate the 
antineutrino bound at the 1 — 3% level. 
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For stronger levels of noise (P=0.8, Figs.(||)) the same comments can be said. The 
structure and position of the allowed regions from combined total rates are practically 
unmodified but the antineutrino yield impose strong restrictions. For an antineutrino 
probability smaller than 3% only some residual, 90% C.L., allowed regions exist at very 
small mixing angle, Am? ~ 10 -6 — 10 -7 eV 2 and moderately high regular magnetic field 
(Plot (C)). The same regions are still acceptable for P=0.7, Figs.(||)). For this level of 
noise something unexpected happens at extremely high regular magnetic field (Plot (D)) 
(probably too high to be acceptable on astrophysical grounds): a new, large, acceptable 
region appear for Am 2 « 1(T 5 . This region disappears again for extremely high chaotic 
fields (P=0.55, Figs.([To|)). This apparently erratic behavior could be worthy of a more 
detailed further study. Note that even for this value of P some residual regions with a 
90% C.L. are marginally acceptable from reconciliation of all experiment total rates and 
antineutrino bounds ||[ if the regular magnetic field is ~ 200 kG ( for fi = 10 _11 //s). 



6 Conclusions. 

We have presented calculations of neutrino spin flavour precession in presence of matter 
and magnetic field for a theoretically motivated solar magnetic field profile. This magnetic 
field, which is the solution to static magnetic hydrodynamic equations, is a twisting field. 
The degree of twisting and the transversal profile are functionally related under this model. 
In our calculations we have considered only the solution with minimum twisting. 

Additionally the effect of a layer of magnetic noise at the bottom of the convective 
zone has been considered, we have justified that the level of noise in this region can be 
certainly very high. 

We have presented expected signals and expected production of antineutrinos with 



a without presence of noise. We confirm previous results [17, |25| for small mixing and 
ad-hoc regular magnetic field profiles. 

We find that MSW regions (Am 2 ~ 10 -5 eV 2 , both small and large mixing solutions) 
are stable up to very large levels of noise (P=0.7-0.8) but they are acceptable from the 
point of view of antineutrino production only for moderate levels of noise (P ~ 0.95). 

This is in agreement with recent results obtained through a direct ensemble averaging 
of the solution to Eq.(|l]) |jl2|j . The stronger r.m.s field occurs at the convective zone, the 
wider (5m 2 , sin 2 29) region should be excluded when considering the constrain imposed by 
existing antineutrino bounds. 

For strong noise, P = 0.7 or bigger and reasonable regular magnetic field, any parame- 
ter region (Am 2 , sin 2 26) is excluded. This model of noisy magnetic field is not compatible 
with particle physics solutions to the SNP. One is allowed then to reverse the problem and 
to put limits on r.m.s field strength, correlation length and transition magnetic moments 
by demanding a solution to the SNP under this scenario. 
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Experiment 


Data ±(stat.) ±(syst.) 


Theory 


SData/SsSM 


Homestake 
Kamiokande 
SAGE 
GALLEX 
SK (504 days) 


2.56 ±0.16 ±0.14 
2.80 ±0.19 ±0.33 
66.61^ 
77.5 ± 6.2it? 

9 ,,+0.05+0.07 
z - 44 -0.09-0.06 


7.7lt;o SNU 
5.15^;? 10 6 cm-V 1 
129±g SNU 
129^ SNU 
5.15±?;?| 10 6 cm-^- 1 


0.33 ± 0.029 
0.54 ±0.07 
0.52 ±0.06 
0.60 ± 0.06 

0.474 ± 0.020 



Table 1: Neutrino event rates measured by solar neutrino experiments, and corresponding 



predictions from the SSM (see Ref.[24] and references therein, we take the INT normal- 
ization for the SSM data). The quoted errors are at la. 



p 




T), i»B = 1) 


T), i»B = 5) 


7?, 0x5 = 20) 


1.00 


0.0 


0.0 


0.00 


0.00 


0.999 


2.2 


5.1 


0.2 


0.01 


0.95 


10 


100 


4 


0.2 


0.80 


23 


540 


22 


1.3 


0.70 


34 


1100 


46 


3.0 


0.55 


51 


2600 


100 


6.6 


0.51 


68 


4600 


180 


11 



Table 2: Values for the noise parameters y (Bq) and r/ (Eq.(||)) assuming Lq = 1000 km. 
All /iB are given in nuB^ units. The quantities P, rj are dimensionless. 
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Figure 1: Transversal magnetic field profile. Solutions Bq± = Bosm9f(r)/r (see Eq. (|Tl|) ) 
with k = 1 (Continuos line) 2, 3, 30 respectively. 
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Figure 2: Probabilities at Earth for a neutrino created at the solar center as a function 
of the magnetic field taken at the bottom of the convective zone (B(0.7ro)), for a fixed 
mixing angle (s| = 0.01). Respectively A) P UpUe , B) P VeV)l , C) P UeW , T))P Ve ^. 
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Figure 4: Relative signal in absence of noise (P = 1) at the different experiments: (A) 
Homestake, (B) Kamiokande, (C) Gallium, as a function of the mixing angle (logs!, hor- 
izontal scale) and squared mass differences (log Am 2 (eV 2 ), vertical scale). The magnetic 
field is fixed to be /j,B = 15 HnB^ at the bottom of the convective zone (r = 0.7 R Q ) 
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;ure 5: The same as Fig.([|) but in presence of noise (P = 0.7, see Table 
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(A) B=1.P=1.00. 
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(B) B=5.P=1.00. 
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(C) B=8.P=1.00. 
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(D) B=20.P=1.00. 
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Figure 6: The parameter regions consistent with the Homestake, Kamiokande and com- 
bined Gallium experiments in absence of noise: P = 1.00. (C.L.= 90%, 95%, 99% (from 
darkest to lighter shaded areas)). Plots (A,B,C,D) correspond respectively to values of the 
regular magnetic field at the bottom of the convective zone: \xB = 1,5,8,20 //.11-B4. The 
electron antineutrino averaged probability, Eq.(|30|), is represented by the dashed lines. 
Present Kamiokande bounds impose the additional restriction P Ve j^ < 0.05 to the total 
rate allowed regions. 
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(A) B=1.P=0.95. 



(B) B=5.P=0.95. 
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(C) B=8.P=0.95. 
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(D) B=20.P=0.95. 
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Figure 7: The parameter regions consistent with the Homestake, Kamiokande and com- 
bined Gallium experiments as Fig.(^) for P = 0.95. 
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(C) B=8.P=0.8. 
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Figure 8: The parameter regions consistent with the Homestake, Kamiokande and com- 
bined Gallium experiments as Fig.iJ]). for P = 0.80. 
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(A) B=1.P=0.7. 
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(C) B=8.P=0.7. 
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Figure 9: The parameter regions consistent with the Homestake, Kamiokande and com- 
bined Gallium experiments as Fig.(^). for P = 0.70. 
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(A) B=1.P=0.55. 
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Figure 10: x 2 contour plots. C.L.= 90%, 95%, 99% (from darkest to lighter), for P = 0.55. 
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